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Abstract 

We apply Dirac's square root idea to constraints for embedded 4- 
geometries swept by a 3-dimensional membrane. The resulting Dirac- 
like equation is then analyzed for general coordinates as well as for 
the case of a Friedmann- Robertson-Walker metric for spatially closed 
geometries. The problem of the singularity formation at quantum 
level is addressed. 

1 Motivation 

One of the most important problems in modern cosmology is the fact 
that according to the Friedmann model the Universe started with the 
singularity. It is expected that once the consistent description of Gen- 
eral Relativity with Quantum Mechanics is obtained the singularity 
will be removed. 

* e-mail: maciej .trzetrzelewski@gmail.com 
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In the 60' DeWitt pQ addressed this problem by assuming that 
the wave functional satisfying the Wheeler-DeWitt equation should 
be zero whenever the singularity occurs at classical level. In the case 
of a Friedmann Universe this implies that the wave function should 
vanish for the scale factor a = 0. DeWitt 's assumption is an example 
of one of several boundary conditions (e.g. the no-boundary condition 
[2], the tunneling condition [3] the symmetric condition [3] - see 
[6] for a review on this subject) which one could set to argue that 
the singularity can be avoided. More recently, it has been pointed 
out that this aim can be also achieved by using techniques of loop 
quantum gravity [TJ, by considering mini-superspace models in D = 11 
supergravity jH] or in (string theory inspired) brane models where one 
addresses this problem from a different perspective in which the Big 
Bang is identified with brane collisions (for a review see [§l I10|). 

A similar singularity formation also takes place in much simpler 
systems of extended objects forming minimal manifolds embedded in 
Minkowski spacetime. These objects will in most cases collapse to a 
point after a finite time (there are exceptions e.g. when they rotate 
which can prevent the collapse: a rotating closed string on a plane 

I12j or a rotating flat torus in M 4 [13 J. In this paper we firstly dis- 
cuss minimal manifolds (Section 2) and then proceed to more general 
geometries (Section 3). 

The idea that one can consider spacetime as a 4-manifold Ai em- 
bedded in some ambient Minkowski space E, and therefore to quantize 
not the metric but the coordinates with give rise to the metric, dates 
back to late 70' |14^15j - see |16j for a comprehensive list of references 
about the subject. Solutions of Einstein's equations give automati- 
cally the solutions for equations of motion for embedded coordinates 
- however the converse is not true [15J. More recently it has been 
argued that the equivalence can be obtained by modifying the embed- 
ding approach (TTJ [181 EH] m a certain way. 

Although it is not clear what is the physical meaning of the em- 
bedding space (i.e. it is not known if there are extra dimensions) let 
us note that there are several advantages of this approach. The most 
crucial used in this paper are related to a clear understanding of 

• the time variable (time on Ai is inherited from E) and therefore 
the Hamiltonian formalism 

• the notion of the probability density - events, such as a 3-sphere 
contraction, take place in E. 
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Here, we would like to develop this idea in a novel context: the 
method used here is an appropriate generalization of our recent work 
|2U| (see also earlier attempts in this direction |2H [2"2~] ) where we have 
undertaken Dirac's idea |23}l24j to formulate a theory of extended ob- 
jects with spin. This approach (different from considering supersym- 
metric version of the bosonic theory) is based on the Dirac's square 
root procedure applied to the constraints for extended objects. By 
writing the Dirac equation for embedded spacetime, in this way, one 
necessarily assumes the existence of the spinor field (actually a func- 
tional) of the manifold considered. Therefore the wave function of 
the embedded Universe will not be a scalar but a spinor. Intuitively, 
this change should make a major difference - after all bosons "like" to 
be in the same state while fermions quite the contrary. Therefore if 
one considers a contracting fermionic 3-sphere, one should observe a 
repulsion effect when the radius of a sphere is sufficiently small. 

A behaviour of this type we find in the case of FRW metric for the 
action containing the Einstein-Hilbert term and the cosmological con- 
stant term. Equations of motion for embedding coordinates are more 
general then Einstein's equation for the metric - in particular they 
allow for singular solutions (at a = 0) in situations where Einstein's 
equations develop no singularities. We find that at quantum level, 
when considering the Dirac equation, these singularities are avoided 
by showing the the wave function is zero for a = 0. 

2 Minimal 4-manifolds 

Embedded 4-manifold A4 will be described here by coordinates X (x), 
A = 0, . . . , D — 1 where x M , [i = 0, 1, 2, 3, is the internal parametriza- 
tion of the manifold. The metric of the manifold is induced from E 
by 

g„u{x) = VAB d^X A (x)d u X B {x) (1) 

where t]ab is a Minkowski metric (throughout the whole paper we are 
using c = h = 1 units as well as the conventions of the Landau-Lifshitz 
textbook [25], in particular the signature of g^ u is (+1,-1,-1,-1) 
and correspondingly r/oo = !)• Such embedding can be found globally 
for many important geometries which are the solutions of Einsteins 
equations |26| . Moreover, it is well known that locally one can always 
find such embedding in 10 dimensional Minkowski space |27j . 

For the action describing the dynamics of X A we take the Dirac 
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action |23t ¥M\ 128] which, for three dimensional membranes, takes the 
form of the cosmological constant term 



S= £\d x, £\ = -Xy/^g, g:=detg^ u . 
Jm 

Here A is a tension and is positive - we will discuss negative tension 
later on. Varying S with respect to X A gives 

{V=99^d v X A ) = 0. (2) 

In addition to the equations of motion, there are constraints 

V A V A = -X 2 detg rs , r,s = 1,2, 3 (3) 

and V A d r X A = satisfied by the canonical momenta 

V A •= dC x /d(d X A ) = -X^gg^d^X A . (4) 

2.1 Dirac equation 

Equation ^ is the counterpart of the mass shell constraint for point- 
like particles. Therefore its linear (in Va) form will serve as a Dirac- 
like equation for extended objects. The square root of the l.h.s. of ^ 
can be easily performed with use of the Dirac gamma matrices ^ A in 
D dimensions. As for the r.h.s. there exist at least two ways of doing 
that. Using the identity 

detg rs = ^{X A ,X B ,X C }{X A ,X B ,X C } 
and a related one 

M 2 = detg TS l, M = ^ ABC {X A , X B , X c } 

where ■y ABC = \{^ AB ^ C + cycl.), j AB = ^b^T" 8 ] and where we used 
the Nambu bracket {f,g,h} = e rst d r fd s gdth, [29], we find that the 
constraint ^ can be linearized as 



[~ ijA 5X^ + \l-li X2 {X A ,X B ,XC}{X A ,X B ,Xc}J * = (5) 

or 

+ ^iX 7 ABC {X A ,X B ,X c }] * = (6) 
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where we substituted the functional derivative Va = —iS/SX . Let us 
note that there is no ordering ambiguity related to functional deriva- 
tives which would have taken place if the equation was second order 
in Va (this is the case in e.g. the Wheeler-DeWitt equation [5j). To 
proceed further a quantum counterpart of the Nambu bracket is re- 
quired. In the case of membranes where instead of a Nambu bracket 
there is a Poisson bracket, one uses the matrix regularization |30J and 
the Poisson bracket is then replaced by —i times the commutator of 
SU(N) matrices. Here one faces the unsolved problem of regularizing 
the volumes. 

Equation ^ seems more elegant due to the lack of the square root. 
Moreover, let us observe that it can be solved exactly by 

* = ^ S = Wb lABC ° J X A{XB,X c ,X D }d*x 

where ^>q is a constant spinor, b is such that r y ABCD ^D = b^ ABC . 
This solution is an analog of similar solutions appearing for mem- 
branes [201 ED E21 E2] . It turns out however that equation ^ gives 
(more) interesting results - a least when one concentrates on spheri- 
cally symmetric case which will be now discussed. 

2.2 Spherically symmetric motion 

Let us now concentrate on the Friedmann-Robertson- Walker line ele- 
ment of spatially closed geometry 

g A Bdx A dx B = dt 2 - a{t) 2 dQ 2 , (7) 

dQ 2 = d X 2 + sin 2 X (d6 2 + sin 2 Odtp 2 ), 
X£ [O,7r],0€ [0,7t] j¥ >€ [0,2tt). 
The corresponding expressions for X A for this geometry are |26j 

X^j^Tut, X*=„( ( >sin X co S e, (8) 

X s = a(t) sin x sin 6 cos <p, X 4 = a(t) sin x sin 8 sin <p 

so that D = 5 - it is a minimal embedding. In order to find a(t) let us 
substitute X A to the equations of motion ^ . We obtain 

3 + 3a 2 + aa = => -at" (a) + 3t'(a) + 3t'(a) 3 = (9) 
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where t = t(a) is the inverse function of a = a(t). Integrating the last 
equation we obtain 



t(a) = t + - 2F1 



1 2 5 

2' 3' 3' 



A 2 a 4 



A 2 



1 



a(t ) 6 (l + a(t ) 2 ) 



where 2-P1 is the Gauss's hypergeometric function. An example of the 
solution for the initial conditions a(0) = 1, d(0) = 1 is given in Figure 
1. We observe a typical behavior, just like in the case of membranes, 



a(t) 




Figure 1: Solution of equation (J9]) with a(0) = 1 and d(0) = 1. 
that a 3-sphere will collapse to a point after a finite time. 



2.2.1 Dirac equation 

Let us now proceed to the quantum theory by linearizig the con- 
straints. Because of the spherical symmetry there is only one dynam- 
ical variable o(t). Therefore the quantum theory will be of quantum- 
mechanical type i.e. we are looking for the equation involving deriva- 
tives w.r.t. a - not the functional derivatives as in Q or @. This 
is completely analogous to considering mini-superspace models intro- 
duced by DeWitt [1] where instead of a functional Wheeler-DeWitt 
equation one considers a partial differential equation involving the 
scale factor. Because of that our starting point are not equations 
^ but the constraints ([3]). 

In the spherically symmetric case the determinant det g rs is simply 
—a 6 sin 4 x sin 2 9 hence the square root of ^ can be written as 

-1 A Va = -XBa 3 sin 2 x sin 9 

where B is a matrix s. t. B 2 = 1 (we choose B = 1), the minus sign 
is a convention. Introducing the canonical momenta pa for the whole 
3-sphere 

PA ■■= / V A dxd9dip 

J t=const. 
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and going to the quantum theory via pa = —id a, we find that 

{-il A d A \a + 2ir 2 \a 3 )i(j(a, t) = (10) 

where | a means taking the a-dependent part of the operator. Therefore 
we have obtained a Dirac equation in (4+l)-dimensions with the cubic 
scalar potential. 



Four comments are now in order. First, equation (10) is a mini- 
superspace counterpart of ([5]), not Q. The mini-superspace equation 
corresponding to ^ would involve the integral of M which turns 
out to be for coordinates ([8]). Therefore one obtains a continuous 
spectrum covering the whole real line in this case. 

Second, note that the choice of the embedding is not unique since 
one could take higher dimensional Minkowski space. However one can 
always choose coordinates for which the momenta Va = for A > 4. 
For such coordinates the index A infl runs from to 4 regardless of 



the embedding hence one arrives at ( 10 ) for any Minkowski embedding. 



Of course other types of geometries in general cannot be embedded in 
(4 + 1) Minkowski space so that D = 5 is not special for the generic 
case, however every Lorentzian 4-dimensional manifold can be locally 
embedded in 10 dimensional Minkowski space [27] - a result which can 
be useful for generalizations of the analysis presented below. 

Third, there is an ambiguity in interpreting \ip(a, t)\ 2 . As a prob- 
ability density of 

• a 3-sphere with the radius £ [a, a + da] at time t 

p(a,t) = \iP(a,t)\ 2 (11) 

• an infinitesimal element d?x of a 3-sphere being at distance G 
[a, a + da] from the origin of the coordinate system at time t- 
therefore the probability density for the whole 3-sphere would 
be 

p(a,t) = ira 3 |VM)| 2 . (12) 
It is in our opinion disputable which one should be chosen although 



the later seems more sensible since equation (10) is written in the 
spherical coordinates hence one expects the measure factor to apear 
(as it is the case in the Dirac or Schrbdinger equation). We will discuss 



both (11) and (12) in this paper 



Forth, from (10) one can consider the evolution of wave packets 



tfj(a,t), by using the Hamiltonian formulation 
idti>{a,t) = Hi(>(a,t) t H = 7 °{ 7 k d k \a + 27r 2 Xa 3 ), k = 1,2,3,4. 
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A detailed analysis of that dynamics will not be discussed in this 
paper. 

The radial part of the Dirac operator in 4 + 1 dimensions can be 
written in general as [M] 



8 a + I pjG = (E+m-V)F, (-d a + ^-j F = (E-m-V)G (13) 

where E is the energy, m is the mass term, V is the potential (the 
Aq component of the gauge field), K = ±(Z + 3/2), I is the angular 
momentum (we shall consider I = from now on). The functions F 
and G are related to the wave function ip by 

V(a, t) = a- 3 / 2 e~ lEt [F(a)Mx, 0, V>) + G(a)Mx, 0, <p)] (14) 

where (pi and (p 2 are suitably chosen, orthogonal, a independent, 
spinors (for more details see |34J). We will normalize them so that 

\^a,t)\ 2 = {F{a) 2 + G{a) 2 )/^a\ 

It is now clear that the choice of the interpretation of \ip(a, t)\ 2 is 
crucial for the analysis of the probability density at a = 0. 

In our case we have m = 2-K 2 Xa 3 , V = 0, therefore introducing 
dimensionless variables x = a(27r 2 | Aj) 1 / 4 , e = £'/(27r2|A|) 1//4 and the 



tension signature a = A/|A| = ±1, the spectral problem (13) can be 
written in the matrix form 

7 / , , ._( -ex 3 d x + f A , . 

V -Ox + f <rx J 

If 4> T = (F,G) is the eigen vector of h m i n fox K = 3/2, A>0 and 
some eigenvalue e m i n then 4> T = (G, F) solves h m i n (fi = —eminfy for 
K = —3/2, A > and 4> T = (—F,G) solves h m i n <p = —e m i n (j) for 
K = 3/2, A<0. Therefore we choose K = 3/2 and A > from now 
on. 



2.2.2 Eigenfunctions at the origin 

Let us start with the following ansatz for functions F and G 

F(x) = x Si e- x / 2 (f + fix + f 2 x 2 + ...), fo^O (16) 
G(x) = x S2 e- x/2 {g + gix + g 2 x 2 + ...), go + 
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i.e. just like in the Dirac equation with the Coulomb potential (in 
which case one obtains s\ = 52)- Depending on the choice of the 
definition of the probability density one obtains different constraints 



on F and G - either F 2 + G 2 (cp. (12)) or (F 2 + G 2 )/x 6 (cp. (11)) 



must be normalizable. Clearly the second choice is more restrictive 
(F and G have to be at least such that (F 2 + G 2 )/x 3 is finite at x = 0) 
therefore to keep the discussion as general as possible we shall assume 
the first constraint. A class of functions obeying the first constraint 
should a priori be larger then for the second one. We shall now prove 



that a careful analysis of the spectral problem (15) actually implies 
that these two classes are in fact equal. 

First, let us use the fact that F and G can be expanded in the 
following orthonormal basis 

«4 a) (*) = -^=x a L^l(x)e-^ 2 , (17) 
T(n + a) 

N n ,a = -, j T\T: Q = mm Si, S 2 , H 6 N+. 

(n + lj! 

where (x) are generalized Laguerre polynomials. Clearly the con- 
dition f^°(F 2 +G 2 )dx < 00 implies that a > however due to the K/x 
term in the operator h m i n one in fact needs to assume that a > 1/2 
- otherwise the matrix representation of h m in would not exist (i.e. 
the scalar products (e n , -e OT ) = J* °° ^e n e m dx would be infinite). This 



constraint can be further improved by using the expansion (16). Be 



cause of the specific form of the operator h m i n (all powers of x are 



integers) we may assume that s\ — S2 G Z. Substituting (16) to (15) 
and analyzing all the possibilities si = S2, si = S2 ± 1, s% = S2 ± 2 
et.c. one finds that there is only one possible choice where a > and 
it is for S2 = si + 1 with si = K{= 3/2). Therefore one should choose 
a = 3/2 which implies that not only F 2 + G 2 but also (F 2 + G 2 )/x 3 
is normalizable. 

Let us note again that the above conclusion is due to the term K/x 
- both the singular behaviour and the coefficient K are important. The 
term comes from the coordinate transformation of the kinetic part of 
the Dirac operator in ( |10[ ) and is unremovable i.e. its generic form is 
K/x with K = ±(2Z + 3)/2 hence it persists for all angular momenta I. 
Analogous situation would not exist for the bosonic field. There the 
corresponding kinetic term would be the Laplace operator which in 
spherical coordinates has the centrifugal barrier term 1(1 — d + 2)/x 2 , 
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in d spatial dimensions. Therefore in the bosonic case there exist two 
sectors of the Hilbert space (for I = and for I = d — 2) where there 
is no singularity. 

We have just shown that \ip\ 2 = (F(a) 2 + G(a) 2 )/a 3 is regular 
therefore if one uses the definition ( |12[ ) then it follows that a 3-sphere 
cannot have a zero radius since the likelihood of that event is zero 
(p(a = 0,to) = 0). On the other hand our analysis also shows that 
\ip(a, t)\ 2 is finite for a = 0. Therefore a similar conclusion could not be 



obtained when taking ( 11 ) as the definition of the probability density. 

Let us finally observe a close analogy between those results and 
the case of Dirac equation in the Coulomb field where \ip(a,t)\ 2 is in 
fact singular while F(a) 2 + G(a) 2 is zero for a = 0. 

2.2.3 The spectrum 



The Hilbert space of the problem (15) consists of square-integrable 
vectors 4> T = (F, G) on [0, oo) satisfying a condition 0(0) = 0. There- 
fore the operator d x is antihermitian in this Hilbert space hence h m in 
is hermitian which proves that the spectrum of h m i n is real. 

Let us also observe that the spectrum of h min is discrete. To see 
this we use an identity 



hence we have an inequality h min > Q . However Q is discrete since 

12 jl _ „2, 



i 2 min > Q 2 - However Q 2 
introducing a± := F ±G the eigen equation Q 2 <f> = r] 2 (J) gives 

h±a± = i] 2 a±, h± := —d 2 + x e ± 3x 2 , 

i.e. a+ and a_ decouple and since h+ and h- are discrete, Q 2 must 
also be discrete. Because Y? min is bounded from below by a discrete 
operator it follows that }^ min itself must be discrete and therefore 
h m in is discrete as well. 

To find the exact spectrum of h m i n we use numerical methods. 
First, we calculate the representation of each operator in h m i n , in 



some orthonormal basis. Our choice is (17) for a = 3/2. Second we 
truncate that infinite matrix representation (i.e. truncate each matrix 
representation of operators appearing in the entries of the 2x2 matrix 
hmin) and then numerically diagonalize the resulting finite matrix. 
The spectra and the eigen vectors of the truncated matrices converge 
to their exact counterparts when the size of the matrix is increased. 
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Figure 2: Convergence of energy levels of h m i n (for K = 3/2, a — 1) with n 
- the size of truncated matrix representation - up to n = 300. 

The results of this numerical approach are presented in Figure 2. It 
follows that the first positive and negative energy levels are 

e«3.4, 5.0, 6.5, 7.8, 9.0, 10.3, 11.4,..., 

-e«2.2, 4.1, 5.7, 7.1, 8.4, 9.7, 10.9,... . 

The plots for F 2 (x) + G 2 (x) are presented in Figure 3. The negative 
energy solutions are normalizable therefore one cannot ignore them. 
This would be troublesome if the 3-sphere was interacting with some 
external field or with another membrane in the embedding space E 
- in which case a 3-sphere could loose energy. Because the spectrum 
is not bounded from below this would imply that the energy of a 
3-sphere could become arbitrary negative. Consistency would require 
that there is no such external fields and that there is only one 3-sphere 
in E. Otherwise one needs to introduce a Dirac sea (of negative energy 
3-spheres) i.e. assume that all the negative energy states are filled. 

3 General action 

In this section we will generalize previous considerations to the case 
of Einstein-Hilbert action with the cosmological constant term and 
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Figure 3: F 2 (x) + G 2 (x) for first three positive and negative energy states of 



hrain- 



matter. Therefore we consider 



S = S\ + Seh + S„ 



(19) 



Sx 



-A 



-g d A x 



Seh 



M 



-gR d A x, 



1 

2k jm 

where k = 8ttG, R is the Ricci scalar, S m is the action for matter. The 
tension A is related to the cosmological constant A via A = A/8irG, 
from now on we will set G = 1. Our view on the action ( |19| ) is 
non-standard here. Usually one takes the Einstein-Hilbert term as a 
starting point and only then adds "something else" e.g. S\ or S m . 
Here we consider space-time as embedded manifold with tension A 
therefore S\ is a starting point while Seh and S m are in addition to 
it. Moreover the Einstein-Hilbert term is defined on A4 not in the 
embedding space E. As such S can be also viewed as an action for 
the vector field X A (x) onM - with no reference to the embedding. 
Varying S with respect to the X A one obtains 



6xS 



M 



-gd x 



where G^ v = R^ v — \g^ y R is the Einstein tensor, R^ in the Ricci 

= 25S m /5g^ u is the energy momentum 



' [ii/ 
tensor (R 



R a 



/ictv. 



2 

and T, 
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tensor. The variation of g^ v gives 2d jJL X^d u X A hence the equations 
of motion are 



g ( -Xg^ + -GT - T» v ) d v X 



0. 



(20) 



These equations can be obtained form the usual Euler-Lagrange equa- 
tions [361137]. 



3.1 Spherically symmetric motion 

If g^u satisfies Einstein's equations with the cosmological constant and 
matter terms 

- XgT + -G^ u = (21) 

K 



then clearly equation (20) is satisfied however it is not so clear if 
there exist coordinates X which give rise to that metric via ([!]). The 
problem lies within the energy-momentum tensor term - if T^ u = 
then one can always find (locally) the embedding coordinates however 
if ^ such coordinates may not exist. We shall now verify this 
remark for the Friedmann- Robertson- Walker line element Q and the 
energy-momentum tensor for a perfect fluid 

T^ u = {p + p)S^qS u0 ~ V9iiv (22) 



There are a priori five equations in (20) for each A however for A 



1,2,3,4 one obtains the same equation which is 

N := 3(l + d 2 ) 2 -3Ka 2 (A-p+(A + p)d 2 )-K(A + p)a 3 a + 9a(l + a 2 )a = 

(23) 

while for A = we find 

N - 3na 2 (p + p) = 0. 

It follows that nontrivial solutions for a exist only if we put p + p = 
hence it is impossible to find embedding coordinates for for FRW 
metric when p > and p > 0. On the other hand if we allow the 
pressure to be negative then the condition p + p = can be obtained 
by considering the matter term S m for the scalar field eft in a particular 
potential V{4>). Then one can show that the equation of state p+p = 
is approximately satisfied for slowly rolling fields (i.e. h<p « V{cj>)) 
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[3H1 [Ml HO] • Let us therefore assume that p = —p. Equation (23) can 
be written in a convenient form 



3 + 3a + ad 



24(1 + a 



2\2 



9 - k(A + /o)a 2 + 9d 2 



(24) 



and we see that in the limit k — > oo one recovers the equation for 
minimal 4- volumes Q. It is now straightforward to verify that (24) 
is a consequence of Firedmann equations for p + p = 0. Substituting 
(22) and Q to (21) we obtain standard equations 



1 , k(\ + p) 



a 
a 



1 



K(\ + p). 



(25) 



Using expressions for a and a obtained from (25 ) we find that the l.h.s. 
and the r.h.s. of (24) coincide (giving 4k(A + p)a 2 /3 each). Therefore 



we have shown that (25) imply (24). It is reasonable to ask if the 



converse is also true i.e. whether (24) are equivalent to (25). We will 



now show that this is not true i.e. the space of solutions of (24) is 



much larger then the space of solutions of ( 25 ) . 

For simplicity let us assume that p = p = i.e. the Universe is 



empty. Then the general solution of the Friedmann equations (25) for 
A > are non singular 



a(t) = 737a K COS h(tv/A/c/3) 



(26) 



(De Sitter solution of spatially closed Universe) while for negative A 
and p = equations ( |25| ) have no solutions. As shown above (26) is 
also the solution of (|24[) however it is not clear if it is a general solution 

of HQ 

Let us first analyze the numerical solutions of (24) for p = - 



Figure 4. We see that not only the solutions do not resemble (26) but 



also that there is a critical value of Xk below which the solutions de- 
velop singularity in finite time (for the boundary conditions a(0) = 1, 
a(0) = 1 the critical value satisfies the bounds 0.056 < Ao^ < 0.057). 
This result is rather interesting since we obtained three classes of solu- 
tions (asymptotically expanding, asymptotically static or contracting) 



1 This is not in conflict with the uniqueness of solutions since the space of boundary 



conditions for (24 1 is much larger then the space of boundary conditions of (25 1. A simple 
example of this sort is a pair of equations a = 1, a = t which have general solution 
a(t) = t. On the other hand they also imply that aa — t the general solution of which is 
a(t) = VWTA. 
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Figure 4: A family of solutions of (24) with a(0) = 1, d(0) = 1 for various 
(positive and negative) values of Xk. There exists a critical value below which 
the solutions develop singularity. 



similar to Friedmann solutions however let us point out that unlike 
for Friedman solutions there is no matter here. Therefore the embed- 



ding equations ( 20 ) allow for solution which do not solve the Einstein 
equations - in particular they allow for singular solutions for = 0. 
Let us now give more rigorous proof that solutions of ( 24 ) are more 
ral 
0) 



general then (26). It is useful to introduce the function F defined as 
(P 



3 d z + l 
k\ a 2 



(27) 



In case of Friedmann equations (25) F is simply 1 - we will now be 
looking for more general F of the form F = F(a). Expressing d in 



terms of F and a and substituting the result to (24) we obtain 



kX (1 - F/3)aF 
1 F-l/3 



(28) 



Equations (27) and (28) are equivalent to (24) for p = 0. Using the 
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formula a = \ t and applying it to ( 28 ) we obtain 



-k\ 



a aF(a)(l-F(q)/3) 
F(a) - 1/3 



On the other hand a from (27) is 



-a z F-l. 



(29) 



(30) 



Equating (29) with (30) we obtain an integral equation for F. The 
derivative w.r.t. a of that equation gives 

1 - F(a) 



aF'(a) = 8F(a) : 



3F(o) - 1 

Integrating the above equation we obtain a non-linear equation for F 

(F - l) 2 F = e A /a 8 , AeR. (31) 

We now see that the Friedmann equations case, F = 1, corresponds 
to taking the limit A — > — oo in (31). However solutions of (31) are 
more general then that. Equation (31) has three roots, two of which 

\ 1/3 



are complex and one is real which is 



F(a) = -(2 + x + l/x), x 



27e A + 3\/3\/27e 2A - 4a 8 e A 



1 



To find the solution for a(t) one can use the definition (27) to arrive 
at the complicated integral 



da 



"0 



'-F(d)d 2 



t - t , a := a(t ). 



This completes the proof that a class of solutions of ( 24 ) is larger then 
a class of solutions of ( 25 ) . 



3.2 Dirac equation 

In this section we shall discuss the singularity formation in the context 
of the Dirac equation using similar method as in Section 2. For general 



coordinates the canonical momenta for the action (19) are now 



n A = v A + (^g 0m - t ^ d^x A 
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where Va are as in Q . The constraint Q for Va still holds therefore 
its linearized form will be 



_ G 0» _ T M \ q Xa 



^ = Ay 7 - det g r 



(32) 



To proceed further we now concentrate on the Friedmann-Robertson- 
Walker line element for closed 3-geometry. Let us therefore assume 
that the metric is given by Q and that there is no matter at all 
(as shown previously this case contains singular solutions). Following 
Section 2 we introduce the average momenta for the whole 3-sphere 



7TA 



HA<ixd9dip 



t=const. 



and, as in (10), we will consider the wave equation by substituting 
ir a — > —id a- Equation (32) becomes now 

{-il A d A \a + 2vr 2 Aa 3 )V(a, t) + j°V^(a, t) 







(33) 



where V is an operator corresponding to the classical expression 
1 f „qo . — . — ^ 2 



V 



G 



-gdxdOdip 



(d 2 a + a) 



where we used the fact that G^ v is diagonal and that the integral 
over d§XA is non zero only if A = (moreover G 00 = 3(d 2 + l)/a 2 



and \f—g = a 3 sin 2 x sin 6). Note that V enters equation (j33h in the 



same way the zero'th component of the electromagnetic field enters 
the Dirac equation with an external field. To find the quantum coun- 
terpart of the term involving the time derivatives of a we first calculate 
the momenta conjugate to a 



Pa 



3vr . 

aa, 

2 ' 



V 



i pI 



3tt 



37r a 



The momenta are calculated using a standard method [5j (in deriving 
p a one uses the action (19) augmented by a boundary term needed 
for the consistent variational procedure - that term does not affect 
the equations of motion nor the constraints). Second, there is an 
ambiguity in choosing the differential operator corresponding to p 2 /a 
term. We choose the ordering following DeWitt pQ as 



Pa 



a l/4^ a l/2^ a l/4- 



17 



Introducing dimensionless units as in previous section we find that 
the eigenvalues of (33) can be obtained from the following spectral 
problem 

hdS<j> = £dS<P, h dS = hmin + v(x)l, (34) 

1 „ 1 „ 1 



v(x) 



~A - ~(3x, 
2/3 2 



x 



l/i dx x l/2 



d x - 



X 



1/4 



where h m i n is as in (15) and /3 = 3^^/ A is a dimensionless constant. 
If we were to substitute the experimental value of A then the term 
proportional to A would be extremely small while the linear term 
would be very large and negative. 



Singularity avoidance 

As in Section 2 we shall now be looking for functions F and G of 



the form (16). The operator A contains terms x~ 3 , x~ 2 d x and x~ l d 2 
which imply that one should take a := min(si, S2) > 3/2 in order to 
have the scalar products (e„, Ae m ) = J °° e n Ae m dx finite. Using the 
same analysis as before we find that this condition can be satisfied this 
time only when si = S2 in which case one finds that the equation for 
(say) si is s\ + si + 7/16 = 0. The roots are 1/4 and 7/4 but a > 3/2 
hence we take a = 7/4 which proves that not only the probability 



density (12) but also (11) is zero for a = 0, concretely 



\^(x,t)\ 2 = {F{xf + G{x) 2 )/^ttx 3 = Vxe- x/2 r(x) 

where r{x) is regular at x = 0. Clearly, this result is much stronger 
then the one obtained in Section 2 - it shows that the eigen functions 
are for a = 0. 

The Hamiltonian formulation is now 

id t ip(a,t) = H%b(a,t), H = ~f°{-f k d k \ a + 2vr 2 Aa 3 ) - V, k = 1,2, 3, 4 

where the wave packets ip(a,tys can be expanded in terms of eigen 
functions of which we know that are for a = 0. Therefore, no 



matter which definition of the probability density one chooses ((11) 



or (12)), one concludes that the 3-sphere cannot have zero radius i.e. 
the likelihood of that event is 0. This certainly fulfills the DeWitt's 
requirement and hence we conclude that the singularity is not present. 

A contracting 3-sphere will approach a minimal size and then start 
expanding. It is likely that this expansion will then be slowed down 
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and stopped due to the x 3 terms in h^s therefore the evolution of the 
wave packet would be cyclic. To verify this remark one would have to 
perform the corresponding numerical analysis. 

The spectrum 

The operator v{x) is negative definite and it enters hds with the plus 
sign therefore the eigenvalues of hds w iU be lowered compared to the 
eigenvalues of h m i n . Moreover has contains a parameter f3 which in- 
fluences the spectrum considerably - that parameter cannot be elim- 
inated by some rescaling of coordinate x. Several spectral lines for 
/3 = 10 and /3 = 20 are presented in Figure 5. 




Figure 5: Spectral lines of has for (3 = 10 (left) and (3 = 20 (right). The cutoff 
reached is n = 1000. To obtain reliable positive spectral values for (5 = 20 
one would have to go to higher n. On the other hand negative eigenvalues 
converge much quicker. 

We observe a discrete spectrum. The positive and negative eigenvalues 
obtained for (5 = 10 are 

era 0.04, 1.2, 2.2, 3.3, 4.3, 5.2,..., 

-era 1.2, 2.6, 6.4, 8.5, 10.4, 12.2,.... 



The eigenvalues, compared to the spectrum of h m i n in ( 15 ), are shifted 
by (roughly) —(3/2. Therefore the lowest positive energy state will now 
have several "bumps" instead of just one bump (see Figure 6). 

The energies are calculated here in units of (2-7r 2 A) 1 / 4 = (ttA/4) 1 / 4 
therefore the density of spectral lines of has and h m in is very large if 
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e=0.04 e=1.2 e=2.2 

F 2 +G 2 F 2 +G 2 F 2 +G 2 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



e=-1.2 e=-2.6 e=-6.4 

F 2 +G 2 F 2 +G 2 F 2 +G 2 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 6: F 2 (x) + G 2 (x) for first three positive and negative energy states of 
h dS for = 10. 

one substitutes /3 ~ 10 60 (i.e. A ~ 1CU 120 ). However the appearance 
of j3 in hds makes the spectrum even more denser (compared to h m i n ) 
so that practically any energy is allowable as if the spectrum was con- 
tinuous. Therefore the smallness of the cosmological constant implies 
in this model that the empty Universe may have arbitrary energy. On 
the other hand if A was say ~ 1 the quantisation of the energies of 
the Universe would be noticeable. 



4 Summary 

By considering spacetime M as as an embedded Lorenzian 4-manifold 
is some ambient, higher-dimensional space E one necessarily faces a 
question whether E is in any sense physical. In this paper we make 
no arbitrary statements concerning this problem however we do take 
advantage of the fact that certain, conceptual problems appearing in 
the formulation of Quantum Gravity by considering Wheeler-DeWitt 
equation (e.g. the problem of time) are naturally solved within this 
approach. Moreover the constraints ^ appearing in such approach 
are very similar, in the form, to the mass-shell constraint for point like 
particles and therefore one is tempted to apply the Dirac's "square- 
root" procedure also in this case. By doing so one arrives at the 
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functional Dirac-like equation for embedded manifolds which are linear 
in momenta. There are at least two distinct ways of doing this 
and (J6|) and in our opinion both are important. 

We then concentrate on the spherically symmetric motions by con- 
sidering 5 3 with one dynamical variable - a(t) - and make a proposal 
(based on Q) for a quantum mechanical system behind such mini- 
superspace model. For minimal 4-volumes we find that the spectrum 
of the corresponding Dirac equation (10) is real and discrete - it con- 
tains positive and negative eigenvalues which are not bounded form 
below. By careful analysis of the wave functions at the origin we 
showed that the amplitude \ip\ 2 is finite at a = hence it is disputable 
whether a 3-sphere can collapse due to the ambiguity of defining the 
probability density. 

It turns out that considering general action with the Einstein- 
Hilbert term one can improve the situation by finding that the ampli- 
tude | tp\ 2 = for a = 0. It is therefore in accordance with DeWitt's 
boundary condition that the function vanishes for those classical con- 
figurations where the singularity appears. This analysis is performed 
for the Friedmann-Robertson- Walker metric for closed spatial geome- 
tries. One would like to call this model an embedded closed Friedmann 
model with no matter (or embedded closed De Sitter model) however 
it is important to note that the classical equations of motion of the 
embedding picture ( 24 ) are much more general then Friedmann equa- 
tions (25). In particular they allow for singular solutions for both 
positive and negative tension A = A/k which is not the case in (25) 
for p = 0. We have shown that in the quantum picture, when one 
considers the Dirac equation (33), these singularities are avoided. 

To generalize this approach one should consider the case p > 
which can be done by introducing slowly rolling scalar fields. 
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